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Abstract 

We determine the specific heat amplitude ratio near a m-axial Lifshitz point 
and show its universal character. Using a recent renormalization group picture 
along with new field-theoretical ei-expansion techniques, we established this 
amplitude ratio at one-loop order. We estimate the numerical value of this 
amplitude ratio for m = 1 and d = 3. The result is in very good agreement 
with its experimental measurement on the magnetic material MnP. It is 
shown that in the limit m — > it trivially reduces to the Ising-like amplitude 
ratio. 
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I. INTRODUCTION 



Among all properties of systems near a second order phase transition, amplitude ra- 
tios of certain thermodynamic potentials above and below the critical temperature play a 
fundamental role together with the critical exponents: they are examples of universal quan- 
tities. All universal amounts are characterized by a limited number of parameters defining 
the universality class and do not depend on the microscopic details of the system under 
consideration. One particular type of universality class is associated to the Lifshitz criti- 
cal behavior [1,2] where the inclusion of competing interactions along one or more spatial 
directions is the main difference with respect to the ordinary critical behavior. 

Recent investigations of the associated critical exponents for the m-axial Lifshitz univer- 
sality class have been put forward using numerical Monte Carlo simulations [3] and field- 
theoretic approaches [4-6]. Within the perturbative e^-expansion, there are two proposals in 
order to unravel the higher loop structure of this sort of critical behavior. One of them makes 
use of a semi-analytic e^-expansion for the critical exponents, where some loop integrals are 
evaluated through numerical integration [4]. Another alternative is the purely analytical 
treatment of all loop integrals involved, such that new renormalization group as well as 6l- 
expansion ideas have been developed in order to determine those critical indices [5,6]. We 
shall focus our attention in the latter thoughout this article for convenience. 

The Lifshitz multicritical point arises in a variety of real physical systems, but we shall 
be concerned here only with its manifestation in magnetic systems where it was originally 
discovered [1]. The uniaxial case m = 1 can be explained in terms of an Ising model with 
ferromagnetic interactions among nearest neighbors spins as well as additional antiferromag- 
netic couplings among the second neighbors along a single axis, known as ANNNI model [7] . 
The competition originates several modulated phases in addition to the usual paramagnetic 
and ferromagnetic phases. It is very simple to analyse the situation near the uniaxial Lifshitz 
point, which arises at the confluence of a modulated and a ferromagnetic phase with the 
paramagnetic phase. Allowing these competing interactions along m spatial directions, one 
obtains the m-axial Lifshitz critical behavior. The anisotropic Lifshitz universality classes 
associated to this sort of critical behavior depend upon the number of components of the or- 
der parameter N, the space dimension d and the number of competing axes m of the system, 
therefore, extending the Ising-like universality classes characterized solely by (N,d). An in- 
teresting feature of the critical exponents of such systems is the property of universality class 
reduction: in the limit m — > the critical indices reduce trivially to those corresponding to 
the usual Ising-like behavior [5,6,8]. 

The theoretical determination of the specific heat amplitude ratio is especially worthwhile 
for two reasons. First, we would like to know whether the universality class reduction also 
holds for this amplitude ratio. As this property is also obeyed by the susceptibility amplitude 
ratio [9] one can conjecture that this property might hold for all critical amplitude ratios. 
Secondly, it would be highly desirable to compare the theoretical value of this amplitude 
ratio with experiments carried out in some real magnetic material exhibiting the Lifshitz 
critical behavior. This would yield a test for the convenience (or not) of the theoretical 
formalism employed in the solution of this problem. From the phenomenological viewpoint, 
theoretical and experimental studies have proved that manganese phosphide (MnP) presents 
a pure uniaxial Lifshitz point (to = 1, d — 3, N — 1) [11, 12]. Moreover, the experimental 



2 



determination of the specific heat amplitude ratio for MnP was realized thirteen years 
ago by means of susceptibility measurements [13]. A theoretical attempt to describe this 
amplitude ratio was performed using the mean field approximation [14] by neglecting the 
contribution of the fluctuations, but the agreement with the experiment was not achieved. 
A proper treatment should include the effect of fluctuations since they play a nontrivial role 
in the determination of this particular amplitude ratio. 

In this paper the specific heat amplitude ratio near an anisotropic m-axial Lifshitz point 
is calculated at one-loop level using a A0 4 field theory combined with new renormalization 
group and e L -expansion methods. We shall treat only the especial Ising-like case N — 1, 
since for iV > 1 and below the Lifshitz critical temperature the appearance of a massless 
Goldstone mode leads to infrared problems which require a separate analysis [15]. In the 
present field-theoretic setting the long-standing difficulty in this calculation is that it requires 
the knowledge of the coupling constant at the fixed point and the specific heat critical 
exponent at two-loop level in order to find out this amplitude ratio beyond the mean field 
approximation. These objects have recently been figured out at two-loop level [6] permitting, 
therefore, the present analysis. It represents the first theoretical report with the effect of 
the fluctuations properly included in the specific heat amplitude ratio for the anisotropic 
Lifshitz critical behavior. We point out that the universality class reduction also holds for 
this amplitude ratio. This is another step forward towards extending this property to all 
other critical amplitude ratios. Furthermore, we estimate it for the uniaxial case m — 1 
in three-dimensional systems. The result is in very good agreement with the experimental 
determination of the specific heat amplitude ratio in the magnetic compound MnP. 

II. SPECIFIC HEAT AMPLITUDE RATIO 

Since the method of calculation has been described in a previous work [9] and is somewhat 
standard for ordinary critical systems [10], we will only review briefly the notations and the 
basic steps. The bare Lagrangian for the m-axial anisotropic Lifshitz critical behavior reads: 

11 111 

L = ^ I Vi 00 I' + V(d-m) 00 I' + S -\ Vm 00 I 2 + 7^000 + 4|Ao0o- (!) 

The Lifshitz critical region is characterized by 5 = with the temperature close but not 
equal to T L . We are going to restrict our analysis using the condition S = henceafter. 

The anisotropic behavior possesses two independent correlation lengths parallel and per- 
pendicular to the competition axes. They allow two independent renormalization group flows 
in momentum space along directions parallel or perpendicular to the competition axes [5]. 
Rigorously speaking, we should assign a label to each renormalized vertex part associated 
to the flow along spatial directions parallel or perpendicular to the competing axes, but 
we have no need to label the renormalized quantities in this work, since the fixed point 
is independent of the flow direction used to renormalize the theory [6]. Consequently, the 
amplitude ratios do not depend on what external momenta scale is varied in order to define 
the corresponding renormalization group transformation. 

The one-loop renormalized Helmholtz free energy density at the fixed point obtained 
by using a nonvanishing quadratic external momenta scale for the m-axial Lifshitz critical 
behavior is given by [16]: 
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F(t, M) = Um 2 + i/M 4 + \(t 2 + g*tM 2 + l -{g*M 2 ) 2 )I SP 

_ i 

where in the above equation t, M (t = Z^~t, <p = Z^ 2 M) are the renormalized (bare) 
reduced temperature and order parameter, respectively, Z^p, , Z^, are normalization functions, 
g* is the renormalized coupling constant at the fixed point, <f is a (d — m) -dimensional wave 
vector perpendicular to the competing axes, whereas A; is a m-dimensional wave vector whose 
components are parallel to the competition axes. The integral I$p is defined by: 



r d d - m gd m k 

ISI ' " J [(((k + K') 2 ) 2 + (q + Py]((k 2 y + q 2 ) ' { } 



We choose to renormalize the theory using normalization conditions where the external 
momenta scale are zero (K' = 0) along the competing axes whereas their nonvanishing 
components are perpendicular to the competition axes. A convenient symmetry point for 
this integral is chosen at P 2 = k\ = 1. As the dimensionful coupling constant is related to 
that dimensionless by g* = (P 2 ) %~u*, where €l = A + ^ — d, this symmetry point transforms 
the dimensionful into the dimensionless coupling constant. The typical geometric angular 
factor for each loop integral characterizing the m-axial Lifshitz behavior is jSd~ m S m T(2 — 
t)T(t)- These factors should be extracted whenever a loop integral is performed and 
absorbed in a redefinition of the coupling constant [5,6,8]. At the symmetry point the 
integral calculated at nonzero external momenta is given by Isp — — (1 + [22]™^), where 

N m = 1 + mi) - V>(2 - f )) and ^(z) = [8]. 

Since the vertex part T^' 2 ^ is additively renormalized, the singular part of the specific 
heat scales with the temperature in the form [5,6]: 

C = -A\t\' aL = -— B(u*) -rg' 2) , (4) 
a L 

where B(u*) is the inhomogeneous term of the renormalization group equation for T^' 2 ^ and 

r (o,2) _ d 2 F(t,M) 

is the vertex part which is related to the specific heat only at zero external momentum 
insertion. 

Above the Lifshitz temperature, M = and using Eqs. (2) and (5) we obtain: 

(0)2 ) _ 1 f d d ~ m qd m k 1 

r ^ ~~2J ( q i + (k^ + ty + 2 Isp - (6) 

When T < T L , we replace the value of M at the coexistence curve, where u*M 2 = —Qt, to 
find 

r (0,2)_ 3 r d d -qd™k 
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The calculation of the one-loop remaining integrals are straightforward. The integration 
over the quadratic loop momenta is trivial, and when the remaining quartic loop integral 
integral is performed we find: 



d d ~ m qd m k t-'-t e L m 



(8) 



where we have used the parameter t in order to unify the language above and below the 
Lifshitz critical temperature T L . 

The exponents v^, cxl and the fixed point (using normalization conditions) were calcu- 
lated in Ref. [6] and expressed as: 



(4 -AO (7V + 2)(iV 2 + 30iV + 56) 2 

aL ~ 2(A7T8) eL ^T7^ ^ 

1 (N + 2 



4(iV + 



4(iV + 8) 3 
1 (N + 2)(N 2 + 23N + QO) 
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u = 



e L { 1 + e L 



(iV + 8) 3 
(9N + 42)" 
(8 + iV) 2 _ 



-Li 



(9) 
(10) 

(11) 



Now we take the particular value N — 1. On the other hand, the additive constant is given 
by: 



B{u*)=K l { 



8T^ 



lp 2 =K 2_ 1 



-^1 £L (1 + h\m£L)- 



After taking k\ = 1 we obtain up to first order in the following result: 



vl2 Df ^ 3 3[i 2 
-B(u ) = h 



a L 



2e L 



1 87 1 n 
H 1 = h2. 

2 4 58 u\ 



(12) 



(13) 



Using Eqs.(6), (7) and (8) in conjunction with the value of I$p we can easily determine r^' 2 -* 
above and below T L . After that, use Eqs.(4) and (12) to derive the specific heat amplitude 
ratio at one-loop level: 



^ = —(i + <d. 



(14) 



where in this formula the specific heat critical index is given at 0(€l, i.e., ai = ^ 



III. DISCUSSION 



Firstly, the result displays a universal form, since it only depends on m and d. Note that 
this easily reduces to the usual critical behavior in the limit m — > 0, therefore confirming 
the reduction of the Lifshitz universality class to the Ising-like universality class. As the 
susceptibility amplitude ratio presents the same property at one-loop level, this might be 
a general feature of all critical amplitude ratios, including not only thermodynamic ampli- 
tudes but also correlation amplitude ratios as well as mixed (correlation/thermodynamic) 
amplitude ratios. Future work will be devoted to those issues. 
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On the other hand, this amplitude ratio was determined experimentally for MnP and 
found to be ^ = 0.65 ± 0.05 [13]. Replacing m — l,d — 3 into the above expression we 

find ^ = 0.74. If we neglect the error bar in the experiment, the difference among the 
two results is about 9%, therefore, in very good agreement with the experimental value 
obtained for MnP. Nevertheless, it is important to mention that this experimental value 
of the specific heat amplitude ratio in MnP was obtained using nonlinear least-square fits 
to adjust it along with the value of the critical exponent oll- The constraint that the value 
of the critical exponents above and below the Lifshitz critical temperature must be equal 
[a\ = consistent with scaling theory) yielded the experimental result oil = 0.46 ± 0.03, 
which in turn led to the value = 0.65 ±0.05 for MnP [13]. Moreover, the different values 
for the exponents produced by unconstrained fits near Tl were recognized by those authors 
as a departure from theoretical predictions, but could also be atributed to the closeness of 
the first-order ferromagnetic- helical (fan) transition, which makes the experimental analysis 
more difficult. Indeed, the experimental specific heat exponent is in disagreement with 
different theoretical estimates for the uniaxial three-dimensional case whose specific heat 
exponent is approximately 0.20 [3-5]. 

Since the experimental specific heat amplitude ratio depends upon the specific heat ex- 
ponent, the best fit which produced the exponent nearly two times larger than that from 
theoretical calculations should affect the value of that amplitude ratio due to error propa- 
gation. Of course, a different fit resulting in another experimental value of does modify 
the above amplitude ratio due to crossover effects which take place when the temperature 
is varied [17]. For instance, when the value of ol is close to 0.270 in a certain fit, the value 
of the corresponding experimental amplitude ratio (fixed by this value of the specific heat 
exponent) is given by £t = 0.435, but the system is at T = 132AK [17] which is away of the 
Lifshitz temperature Tl = (121 ± 1)K for MnP and is affected by the crossover to the Ising 
universality class. The funny surprise is that very close to the Lifshitz critical temperature 
this error is rather small when theoretical and experimental estimates for this amplitude 
ratio are compared as discussed in the present work, in spite of the large deviation between 
the measured specific heat critical exponent for MnP and calculated values using theoretical 
tools as discussed above. 

A comparison of theoretical and experimental results of this critical amplitude ratio with 
the one associated to the Ising-like universality class is illuminating. For N = \,d — 3 the 
e-expansion (e = 4 — d) result of this amplitude ratio at one-loop level yielded = 0.56 [15]. 
By neglecting the error bars, specific heat measures in different materials belonging to the 
same universality class resulted in critical amplitude ratios varying from 0.53 to 0.58 [18]. 
On the other hand, series expansion studies of this ratio produced values in the range 
0.62 — 0.70 [19]. The e-expansion estimate furnishes better values than those coming from 
series expansion even though the e parameter is not a small number in this case (e = 1). 

For three-dimensional uniaxial systems €l = 1.5 is even bigger than in the situation 
for the Ising-like case analysed above. This might be the main consequence for a greater 
difference between the experimental and theoretical results of this amplitude ratio when 
compared to the ordinary critical behavior. In order to see whether this is a good argument, 
we suggest the measurement of this specific heat amplitude ratio in other uniaxial magnetic 
materials like CoNbiO^ [20]. For Ising-like magnets the deviation of these specific heat 
amplitude ratios for different materials is rather small among themselves as well as when 
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they are compared with results derived from the e-expansion. Hence, from our e^-expansion 
results we would expect a larger deviation of this amplitude ratio for different materials 
presenting a uniaxial critical behavior. Nevertheless, our results are rather encouraging in 
comparison with experiments in MnP, indicating the reliability of the e^-expansion approach 
to this sort of behavior. Since further numerical works have not been pursued for critical 
amplitude ratios pertaining to the anisotropic Lifshitz universality classes, we hope our 
endeavor can stimulate other numerical estimates like series expansions or Monte Carlo 
simulations, in order to improve the comprehension of as many critical amplitude ratios as 
possible for competing systems. 

In summary, we have derived the specific heat critical amplitude ratio above and below 
the Lifshitz critical temperature for a general anisotropic behavior 1 < d < m — 1. The 
anisotropic universality class reduces easily to the Ising-like one in the limit m — > 0. Needless 
to say, the result is independent of the choice of the symmetry point employed: had we 
started with a symmetry point characterized by zero external momenta perpendicular to the 
competing axes and nonvanishing external momenta components parallel to the competition 
axes, we would have arrived at the same result. The uniaxial value of this ratio is in very 
good agreement to that measured for MnP. Hence, the analytical approach presented here 
in order to treat this amplitude ratio is convenient to describe phenomenological aspects 
of real physical systems presenting competing interactions. A thorough study of all critical 
amplitude ratios as well as a nonperturbative proof of their universal character based solely 
on renormalization group arguments for both anisotropic and isotropic cases will be an 
interesting topic for future investigation. 

The author thanks C. C. Becerra for useful discussions and financial support from 
FAPESP, grant number 00/06572-6. 
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